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Abstract 



O 

o 

In the framework of quark models, hybrid mesons are either seen as two-body qq systems with 



an excited flux tube connecting the quark to the antiquark or as three-body qqg systems including 
a constituent gluon. In this work we show that, starting from the three-body wave function of 
the qqg hybrid meson in which the gluonic degrees of freedom are averaged, the excited flux tube 



the constituent gluon approach is confirmed for heavy hybrid mesons but, for the first time, it is 



picture emerges as an equivalent qq potential. This equivalence between the excited flux tube and 

O 

shown to hold in the light sector too, provided the contribution of the quark dynamics is correctly 

o ' 

■ taken into account. 
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The study of hybrid mesons is an active domain in theoretical and in experimental particle 
physics. From a theoretical point of view, they are interpreted as mesons in which the color 
field is in an excited state. Numerous lattice QCD calculations have been devoted to hybrid 



mesons 



ELB 



2|, as well as many studies involving effective models. In particular, within the 
framework of quark models, there are two main approaches. In the first one, the quark and 
the antiquark are linked by a string, or flux tube, which is responsible for the confinement. 
In this stringy picture, it is possible for the flux tube to fluctuate, and thus to be in an 
excited state 0, 4]. The second approach assumes that the hybrid meson is a three-body 
system formed of a quark, an antiquark, and a constituent gluon. Two straight strings then 
link the gluon to the quark and to the antiquark. This picture has been firstly studied in 
Refs. [5J, but also in more recent works 

It was suggested in Ref. 8| that, in the static quark limit, the constituent gluon picture 
is equivalent to the excited flux tube one, the total energy of the constituent gluon being 
equal to the energy contained in the excited string. These results were extended in Ref. {9], 
where the dynamics of the quarks has been taken into account. For further developments, 
it is useful to briefly sum up the key points of this last reference. 

Assuming the Casimir scaling hypothesis, it can be shown that the flux tubes in a qqg 



system are two straight strings linking the gluon to the quark and to the antiquark jj], in 
agreement with Refs. [ij]. In this case, taking only the confining interaction into account, 
we can write a spinless Salpeter Hamiltonian for the system, 



i=q,q,9 3=9>Q 

with m g = 0. The three-body eigenequation 

H 3b V 3b {r, y) = M 3b V 3b {r, y) (2) 

can be analytically solved by using the auxiliary field technique if the quark and the antiquark 
are of the same mass. We will only focus on this case in the following. It is worth noting that 
the quark-antiquark separation, and that y is the second relative variable, 
directly linked to the gluon position. Various approximations are necessary to perform this 
resolution, and lead to eigenfunctions which are separable, i.e. 

* ab (r,y)=A(r)B(y). (3) 
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We stress that this exact separability is only an artifact of the auxiliary field technique. 

It can be shown that, by dropping the "gluonic" part B(y), the quark-antiquark wave 
function A{r) = ^2b(f) satisfies the eigenequation 

H 2b ^2b{r) = M 3b ^ 2b (r), (4a) 

where i?2fc is the two-body spinless Salpeter Hamiltonian 



^» = Z)v^ +T ^ +v *^' (4b) 

i=q,q 

with V q q the equivalent two-body potential. 

For heavy quarks, this potential has the form j^] 



V*{r) = v/crV + 27ra(jV + 3/2), (5) 

where a is a new string tension. In this potential, Af = 2n y + £ y , with n y and £ y the 
radial quantum number and the orbital angular momentum with respect to the variable 
y. Consequently, Af defines the gluon state, and the equivalent potential depends on this 
quantum number. Actually, formula (jSJ) is nothing but the energy of an excited string of 
length r [h]], whose square zero point energy is given by 3na (for J\f = 0). It is different 
of the generally accepted value in string theory, which is —2ira(D — 2)/24, with D the 
dimension of space (see for example Ref. [llj, p. 231]). Together with D = 26, it ensures 
that the Lorentz invariance is still present at the quantum level. However, the string we are 
dealing with is an effective one at D = 4, simulating the confining interaction for an excited 
color field. In that sense, our nonstandard value is more relevant for the study of hybrid 
mesons, since 3na is actually equal to the square zero point energy of the gluon and the 
two strings in the qqg system jjj. It is also important to notice that the string tension a is 
not necessarily equal to the string tension a, since the excited string is an effective object 
emerging from the gluon-plus-string system. 

In the limit where the quarks are massless, computations are more complex, and only an 
asymptotic approximate form can be computed for the equivalent two-body potential, that 
is Q 

V l qq -{r > a' 1 ' 2 ) ^ar + -(M + 3/2). (6) 

The purpose of this work is to extend the results of Refs. 8, ^ by performing an accurate 
numerical resolution of the eigenequation ([2]) followed by an accurate numerical inversion of 



the eigenequation (jlj). It is then possible to compute the equivalent potential V qq - for the qq 
pair contained in the qqg system from the corresponding internal qq wave function. These 
results, obtained without the approximations of the auxiliary field method, will be computed 
for hybrid mesons formed of heavy as well as for light quarks of the same flavor, and then 
compared to Eqs. (j5J) and ([6]). The procedure is the following: 

Firstly, the mass M 3b and the three-body wave function Nt^r, y) have to be computed. 
For the lowest lying state, the wave function reads 



*3&(r,!/) 



[3,3] 8 ,8 



lui] 1 ,0 



[1/2, l/2] s ™ , lj ®$£=o(r,i/), (7) 

where the color, isospin, spin, and space functions have been explicitly written. The color 
function is unique since the hybrid meson is in a color singlet, and the isospin function 
is trivial. The spin functions are unambiguously given when S = or 2 since S qq = 1. 
Furthermore, as the C-parity, C = (— j s a good quantum number [5|, there is no 
coupling between 5=1 states with either S qq = or 1. The spin function is thus always 
unique. As we deal with the ground state of Hamiltonian (1T1), only the L = spatial wave 
function has to be computed. It means that the parity P of the state is always positive jlj]. 
For this space wave function, Gaussian trial functions of the following form are used 

JV 

®L=o(r, y) = ^2ai exp (-Oj r 2 - h y 2 - 2 a r ■ y) , (8) 

8=1 

where the variational parameters Oj, fej, and have to be determined by a minimization on 



the mass (see Ref. [l2j for more details on this method). Calculations are simplified 
if we set q = in Eq. (jSJ), which means that the orbital angular momenta relative to the 
Jacobi coordinates r and y are both equal to zero. It can be checked that this has a little 
influence - around 1% - on the ground state mass M^. As Hamiltonian is spin- and 
isosp in-independent, the states with S — 0, 1, 2 and I = 0, 1 (for light quarks) are degenerate. 
Consequently, the ++ , 1 ++ , l +_ , and 2 ++ states are degenerate ground states of our model. 

Secondly, the two-body wave function i &2b{r) has to be computed from ^/ 3 b(f,y) (up to 
an irrelevant normalization factor). Due to good C and P quantum numbers, the color, 
isospin, and spin functions of the qq pair are uniquely determined from ^^(r, y). A general 
way of averaging the space gluonic degrees of freedom is to solve the following equation 

*»(*■) = / f(yM= (r,y)dy, (9) 



where n is an arbitrary exponent, and f(y) an arbitrary weight function. If <E>2 =0 (r, V) were 
separable as in Eq. ([3]), any choice of f(y) and n would give identical results. This is not 
a priori the case for the trial wave function (jHJ). We use in the following the two simplest 
prescriptions, both characterized by f(y) = 1. We call the first one the linear prescription, 
with n — 1 and 

N ( 71 \ 3/2 

d ex p(-^ r2 )' ( io ) 

i=l V «/ 

where r is the quark- ant iquark distance. This prescription is the most straightforward way to 
remove the space gluonic part from the three-body wave function. The second prescription 
we suggest is called the quadratic prescription, for which n = 2, 

N ( 7T \ 3/2 

i,j=l ^ 

It has the interesting feature that J ^^(r) dr — 1 is a priori satisfied when ^3;, has a unit 
normalization. In general, when taking the square root, the sign has to be chosen such that 
^2b{f) is smooth, but no such choice has to be done here since we deal with the ground state 
wave function. The regularized qq wave function is then given by u(r) = r^ 2 b(^)- 

Thirdly, knowing the two-body wave function u(r), the mass M3&, and the quark masses 
(same values as in Hamiltonian H^), the equivalent potential V q q(r) satisfying Eqs. (jlj) can 



)e computed. We gave in Ref. 13| a procedure, relying on the Lagrange mesh method 
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151 ] . to make such computations with a semirelativistic kinematics. This procedure has 



been shown to be accurate up to 3% in the computation of the equivalent potential, and has 
already been applied to study the effective potential between two gluons If]]. 

The procedure we have described can now be applied to a physical system. Let us begin 
with the ccg hybrid, for which we use the parameters of Ref. m c = 1.3 GeV and a = 0.175 
GeV 2 . These values are fitted so that the mass of the J/^> meson is correctly reproduced 
as well as the Regge slope of the p meson family. We begin by showing the evolution of the 
equivalent potential with the number of Gaussian functions in Fig. [U It appears that the 
shape of the potential does not change anymore when N > 6. Moreover for N = 7, the error 
on the ccg mass is less than 1 MeV, which is a very satisfactory precision. Although we only 
plotted the effective potential computed with the linear prescription, we checked that the 
same conclusions hold for the quadratic one but also for other choices. 

The equivalent cc potential is then computed with N = 7, with the linear and quadratic 
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FIG. 1: Evolution of the equivalent cc potential for different values of N. The linear prescription 
and the parameters of Table U are used. 

TABLE I: Numerical values for the different parameters involved in our computations of qqg 
systems, a = 0.175 GeV 2 , m g = 0, m n and m c are taken from Ref. Table 2]. The symbol n is 
used for both the u and d quarks. 



Quark b c s n 

m q (GeV) 4.567 1.300 0.180 0.000 
M 3b (GeV) 10.665 4.323 2.516 2.421 
a (GeV 2 ) 0.200 0.200 0.207 0.207 
W q (GeV 2 ) 0.000 0.160 0.460 0.530 
H q (GeV) 4.757 1.409 0.435 0.374 

prescriptions, and plotted in Fig. [2J It can be seen that both prescriptions roughly lead to 
the same potential. Furthermore, we actually checked that the equivalent potential is nearly 
independent of the prescription which is used. This means that the dynamics of the gluon 
is quasi-decoupled from the quark- ant iquark pair in the ccg system, and more generally 
in qqg systems. The equivalent potential can be fitted by the following generalization of 
potential (jSJ) 

Vqq(r) = y / a 2 r 2 + 2^(^ + 3/2) + W q} (12) 
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FIG. 2: Effective potentials and regularized wave functions (with an arbitrary normalization) 
computed with N = 7 for the ccg system with the quadratic prescription (solid line) and the linear 
one (dotted line) . Potential (fl2|) is plotted for Af = (dashed line) ; it is practically indistinguishable 
from the result of the linear prescription on this figure. The various parameters are summed up in 
Table E 

where Af = 0. It is checked that the introduction of a constant W q ^ in potential (fl2i) is 
necessary to recover the numerical data (see Tabled]). 

Very good descriptions of mesons and baryons, even light ones, are obtained within the 
framework of potential models. It has been recently shown that such models give gg and 
ggg glueball spectra very similar to results predicted by lattice QCD works [171 ] . So we can 
assume that the procedure developed in this paper to study hybrid mesons remains valid for 
light quarks. We thus computed the effective potential V qq for all ground state qqg systems. 
Results are presented in Fig. [HI Let us note that the masses of the s and b quarks were 
fitted on the masses of the 0(1020) and of the Y(IS') by using the model of Ref. [7]. It 
appears from a fit that the equivalent potentials are all compatible with the form (I12p . with 
o m 0.200 ± 0.010 GeV 2 , and a value of W q which increases as the quarks become lighter 
(see Table U). 

As we mentioned before, a 7^ a since the emergent excited string has not to share the 
same string tension than the two initial strings in the qqg system. Such a feature could 
not be found in the previous study of Ref. [9| due to the approximate method used. It has 



been shown in Ref. 18j that the energy of a string in its ground state is affected by the 
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FIG. 3: Effective potentials computed with N = 7 and the linear prescription for the different qqg 
systems (symbols). Potential (fT2|) is plotted for M = (solid lines). The various parameters are 
summed up in Table HI 

presence of dynamical particles at its ends. In the same way, we argue that W q could be a 
modification of the square zero point energy of the two strings because of quark dynamics. 
A first hint is that Wb = is compatible with the data. This corresponds to the static 
limit of potential (jSJ), in agreement with lattice QCD [2]. For the first time however, we are 
able to show that the excited string picture also holds for light quarks, but with a nonzero 
value for W q . A comparison of potential (fT2l with the approximate result © for r ^> a -1 / 2 
leads, for J\f = 0, to W n ~ (4 — 7r)3cr m 0.52 GeV 2 in the limit of massless quarks, in 
agreement with the fitted value. Moreover, defining the dynamical mass of the quark as 
jj, q = (^36 1 \lv\ + m 1 |^3b) 3 j ^ can be checked that W q ~ 0.2//i 9 reproduces very well 
the fitted values given in Tabled! Consequently, W q is directly related to the quark dynamics 
at the end of the strings. 

In conclusion, we have shown that, in the framework of a potential model with spinless 
constituent quarks and gluons, it is possible to compute the dominant interaction between 
a quark and an antiquark in a hybrid meson from a qqg system. The potential obtained is 



a, 



in good agreement with results predicted by lattice QCD calculations for heavy quarks 
provided the string tension is slightly greater in hybrid mesons than in ordinary mesons. 
Moreover, up to our knowledge, we predict, for the first time, that the dominant quark- 
antiquark interaction in a hybrid meson depends on the quark flavor just by a square shift 



S 



in energy, due to the quark dynamics. 

Since our method can easily be generalized to spin-dependent Hamiltonians, we think 
that these results open the way for a deeper understanding of the fine structure of the QCD 
string spectrum: By adding proper short-range interactions [19], the equivalent potentials 
wiH al80 depend on S Qi through spin-orbit or spin-sph, couplings. These eoup.ings could 
explain the short-range splittings of the potentials as they are observed in lattice QCD [2j. 
We leave such studies for future works. 

C. S. and F. B. thank the FNRS for financial support. V. M. thanks the IISN for financial 
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